Two-color photoassociation spectroscopy of ytterbium atoms and the precise 
determinations of s-wave scattering lengths 
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By performing high-resolution two-color photoassociation spectroscopy, we have successfully de- 
termined the binding energies of several of the last bound states of the homonuclear dimers of six 
different isotopes of ytterbium. These spectroscopic data are in excellent agreement with theoretical 
calculations based on a simple model potential, which very precisely predicts the s-wave scattering 
lengths of all 28 pairs of the seven stable isotopes. The s-wave scattering lengths for collision of two 
atoms of the same isotopic species are 13.33(18) nm for 168 Yb, 3.38(11) nm for 170 Yb, —0.15(19) nm 
for 171 Yb, -31.7(3.4) nm for 172 Yb, 10.55(11) nm for 173 Yb, 5.55(8) nm for 174 Yb, and -1.28(23) 
nm for 176 Yb. The coefficient of the lead term of the long-range van der Waals potential of the Yb2 
molecule is C 6 = 1932(30) atomic units {E h al « 9.573 x 10~ 26 J nm 6 ). 

PACS numbers: 34.50.Rk, 34.20.Cf, 32.80.Pj, 34.10,+x 
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I. INTRODUCTION 

A collision between two atoms is a fundamental phys- 
ical process and can be described by a few partial scat- 
tering waves for ultracold atoms. At a sufficiently low 
temperature the kinetic energy of colliding atoms be- 
comes less than the centrifugal barrier, and only the s- 
wave scattering is possible. The s-wave scattering length 
is an essential parameter for describing ultracold colli- 
sions. It also governs the static and dynamic properties 
of quantum degenerate gases like a Bosc-Einstcin conden- 
sate (BEC) or a degenerate Fermi gas (DFG) of fermionic 
atoms in different spin states. Since the s-wave scattering 
length is very sensitive to the ground state interatomic 
potential, especially at short internuclear distance, the 
precise ab initio calculation of the scattering length is 
very difficult, and therefore we must resort to experi- 
mental determination. The most powerful approach for 
determining the scattering length is to measure the bind- 
ing energy (Eb) of the last few bound states in the molec- 
ular ground state, since the energy E\> is closely related 
to the s-wave scattering length. So far, the binding ener- 
gies were measured via two-color photoassociation (PA) 
spectroscopy for Li 0], Na @, K Q, Rb [|, Cs @ and 
He @ . A schematic description is shown in Fig. [IJ If a 
laser field L 2 is resonant to a bound-bound transition, it 
causes an Autler-Townes doublet. This effect is detected 
as a reduction of a rate of a free-bound PA transition 
driven by the other laser field L\. This scheme is called 
Autler-Townes spectroscopy, and it is also explained in 
terms of the formation of a dark state. If both lasers 
are off-resonant and the frequency difference matches E^, 
these lasers drive a stimulated Raman transition from the 
colliding atom pair to a molecular state in the electronic 



ground state. In this Raman spectroscopy, the resonance 
is detected as an atom loss. 

The s-wave scattering length of two colliding atoms is 
determined by the adiabatic Born-Oppcnhcimcr interac- 
tion potential V(r) between the two atoms, which is very 
well approximated at large interatomic separation r by 
the van der Waals contribution, V(r) as —Cq/t 6 , where 
Cq is the van der Waals coefficient due to the dipole- 
dipole interaction. The s-wave scattering length is given 
by the following formula, based on a quantum correction 
to the WKB approximation so as to be accurate for the 
zero-energy {E = 0) limit 0, BQ, 



1 - tan $ - 



Here a 



(1) 



= 2~ 3 / 2 ^|^ (2/xC 6 /^ 2 ) 4 is a characteristic 



length associated with the van der Waals potential, where 
r is the gamma-function, fi is the reduced mass, and % 
is the Planck constant divided by 27r. The semiclassical 
phase $ is defined by 



$ = 



y/-V(r)dr, 



(2) 



where r$ is the inner classical turning point of V(r) at 
zero energy. The number of bound states N in the po- 
tential is [8J 



N = 



(3) 



where [. . .] means the integer part. As was mentioned 
above, the scattering length is very sensitive to the phase 
$ and can take on any value between ±oo when $ varies 
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FIG. 1: (color online) Schematic description of the two-color 
PA spectroscopy. The laser L\ drives one-color PA transition. 
The laser L2 couples the bound state in the excited molecular 
potential to the one in the ground molecular potential. The 
detuning A of the PA laser with respect to the one-color PA 
resonance is set to several MHz for the Raman spectroscopy, 
while A is set to zero for the Autler-Townes spectroscopy. 



over a range spanning n. While the accurate calcula- 
tion of <!> is very difficult since it requires a knowledge of 
the whole potential, <E> is proportional to ^/Ji, and so the 
formula in Eq. JTJ) gives a simple mass-scaling of the iso- 
topic variation of the scattering length once a and N are 
known for one isotopic combination. Mass scaling applies 
as long as small mass-dependent corrections to the Born- 
Oppenheimer potential [III, [III, 03 can be ignored. While 
mass scaling is often used for bound and scattering states 
for cold atomic systems [hJ El, and seems applicable 
to large-mass systems within experimental uncertainties 
fljj ]. exceptions are known [Ttij ]. and its accuracy should 
be carefully tested for different kinds of systems. 

Ytterbium (Yb) is a rare-earth element with an elec- 
tronic structure similar to that of the alkaline-earth 
atoms. One of the unique features of Yb atoms is a 
rich variety of isotopes with five spinless bosons ( 168 Yb, 
170 Yb, 172 Yb, m Yb, 176 Yb) and two fermions ( 171 Yb 
with the nuclear spin 1=1/2 and 173 Yb with I = 5/2), 
which enables us to study various mixtures of dege nerate 
gases of Yb atoms. In fact, recently, a BEC EJ and a 
DFG [IH of Yb atoms have been achieved. Another dis- 
tinct feature of Yb atoms is the simple electronic ground 
state of 1 S'o symmetry. Therefore, the ground molec- 
ular state of Yb has only one potential of 1 S g molec- 
ular symmetry with no electronic orbital and spin an- 
gular momenta. This is in contrast to the case of an 
alkali dimer which has spin-singlet 1 Y> g and spin-triplet 
3 E ff ground states with a complicated hyperfine struc- 
ture. These two unique features of the Yb system, that 
is, the existence of rich variety of isotopes and one sim- 
ple isotope-independent molecular potential, offers the 
possibility to systematically check the scattering length 
theory with an unprecedented precision. If we can be 
confident in the theory, then it would enable us to deter- 
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FIG. 2: (color online) Schematic description of an experi- 
mental setup for PA lasers. The two lasers were prepared by 
splitting one laser, and two double-pass acousto-optic modu- 
lators (AOMs) were employed to tune them. Then the two 
lasers were aligned in the same path. The solid lines and 
broken lines indicate the laser beams before and after the 
double-pass, respectively. A polarizing beam splitter (PBS) 
and a A/2 plate after the optical fiber were employed to fix 
polarization of both lasers in the same direction. 



mine the scattering lengths of all possible isotope pairs 
which are not measured experimentally. So far, to deter- 
mine the s-wave scattering length of 174 Yb, one-color PA 
spectroscopy was performed [l9( , which gave the result of 
5.53(11) nm. From the cross-dimensional rethcrmaliza- 
tion technique, the absolute value of the scattering length 
for 173 Yb was estimated as 6(2) nm [III, assuming that 
the spin was completely unpolarized. However, these two 
values are not enough for a rigorous test of the theory. 

In this paper, we report an accurate determination 
of the s-wave scattering lengths for all Yb isotopes in- 
cluding those for different isotope pairs. By using two- 
color PA spectroscopy with the intercombination transi- 
tion 1 5o— 3 Pi, we successfully determined the binding en- 
ergy Eb of twelve bound states near the dissociation limit 
of four homonuclear dimers comprised of bosonic atoms 
( 170 Yb 2 , 172 Yb 2 , 174 Yb 2 , and 176 Yb 2 ) and two dimers 
comprised of fermionic atoms ( 171 Yb 2 and 173 Yb 2 ). The 
spectroscopically measured binding energies are in excel- 
lent agreement with theoretical calculations based on a 
simple model potential that was fit to the data. The cal- 
culated s-wave scattering lengths for the six isotopes obey 
the mass-scaling law with very good precision. Moreover, 
this excellent agreement allows us to accurately deter- 
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FIG. 3: (color online) Two-color PA spectra at about 1 fiK. The horizontal and vertical axes are the frequency difference 
between the two lasers and the number of remained atoms, respectively. The vertical error bars represent the fluctuations 
of the number of atoms. The solid lines are fits of the Lorentz functions. The Autler-Townes spectroscopy is applied to the 
v = 1, J = 0, 2 states of 171 Yb2 and the v = 1, J = state of 172 Yb2. while the Raman spectroscopy is applied to the others. 



mine the scattering lengths of all twenty-eight different 
isotopic combinations. In addition, we can reveal scat- 
tering properties of other partial waves such as p- and 
d-wavc scatterings and energy dependence of the elastic 
cross sections. These results are an important foundation 
for future research, such as the efficiency of evaporative 
cooling, stability of quantum gases and their mixtures, 
and the clock shift [2fj |. 



II. EXPERIMENT 

The experimental setup was almost the same as our 
previous experiment of one-color PA spectroscopy [2ll ]. 
All the experiments were performed at about 1 /zK, 
where only the s-wave scattering is possible. Atoms 
were first collected in a magneto-optical trap (MOT) 
with the intercombination transition -"-So — 3 Pi at 556 
nm. The lincwidth and saturation intensity of the transi- 
tion were 182 kHz and 0.14 mW/cm 2 , respectively. The 
laser beam for the MOT was generated by a dye laser 
whose linewidth was narrowed to less than 100 kHz. The 
dye laser was stabilized by an ultralow expansion cavity, 
whose frequency drift was typically less than 20 Hz/s. 
The number, density and temperature of atoms in the 
MOT were about 2 x 10 7 , 10 11 cm" 3 and 40 ^K, respec- 
tively. Then the atoms were transferred into a crossed 



far-off resonant trap (FORT). The number, density and 
temperature of atoms in the FORT were about 2 x 10 6 , 
10 13 cm -3 , and 100 fiK, respectively. To reach lower tem- 
peratures, evaporative cooling was carried out by grad- 
ually decreasing the potential depth of the horizontal 
FORT beam to several tens of fj,K in several seconds. 
Evaporative cooling worked rather well for the bosonic 
isotopes, 170 Yb, 172 Yb and 174 Yb and for the fermionic 
isotope 173 Yb. Typically 1 x 10 5 atoms at the temper- 
ature of about 1 /iK finally remained in the trap and 
the density was between 10 13 cm -3 and 10 14 cm -3 , al- 
though optimized evaporation ramps and efficiency for 
each isotope were different. In particular, we observed 
rapid atom decay for 172 Yb due to three-body recom- 
bination. The steeper evaporation ramp was needed for 
172 Yb to obtain enough number of atoms. It is also noted 
that an unpolarized sample of 173 Yb enabled us to per- 
form efficient evaporative cooling even at a low temper- 
ature via elastic collisions between atoms with different 
spins. For the bosonic 176 Yb and fermionic 171 Yb iso- 
topes, however, the evaporative cooling did not work well. 
In order to cool 171 Yb and 176 Yb, we performed sympa- 
thetic cooling with bosonic 174 Yb [13, [23[. Bichromatic 
MOT beams for simultaneous trapping of two isotopes in 
the MOT were generated by an electro-optic modulator 
(EOM), of which the modulation frequency corresponds 
to the isotope shift. Typically 6 x 10 4 atoms for 171 Yb 
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and 2 x 10 5 atoms for 176 Yb at the temperature of about 
1 [iK finally remained in the trap. The bosonic isotope 
of 168 Yb was hard to collect in the MOT due to its small 
natural abundance of 0.13 percent, although there are no 
fundamental difficulties for 168 Yb in principle. 

After the evaporative cooling, the two lasers, L\ for the 
free-bound transition and Li for the bound-bound tran- 
sition, were simultaneously applied to the atoms in the 
trap for about 30 ms. These beams were focused to a 100 
/Ltm diameter. The schematic setup for the PA lasers is 
represented in Fig. [2] The two laser beams were prepared 
by splitting one laser beam for the MOT and therefore 
had the same frequency linewidth and stability as the 
MOT beam. The relative frequency was controlled by 
acousto-optic modulators (AOMs) . The two laser beams 
were coupled in the same optical fiber. A polarizing beam 
splitter (PBS) and a half- wave (A/2) plate after the op- 
tical fiber were inserted to fix the polarization of both 
lasers in the same direction. Finally, the PA laser beams 
were aligned to pass through the atoms in the FORT by 
using a CCD camera for absorption imaging. The detun- 
ing of the PA laser with respect to the atomic resonance 
1 Sq — 3 Pi was easily checked by observing the frequency 
at which the atoms in the MOT disappeared. The power 
of the PA laser was monitored by a photodiode. For the 
Raman spectroscopy, the frequency fi of L\ was fixed 
with a certain detuning A from a particular PA reso- 
nance, and the frequency fi of Li was scanned to search 
for the bound states of the ground state. For the Autler- 
Townes spectroscopy, f\ was fixed to a particular PA res- 
onance, and fi was scanned. Table U shows the excited 
state rovibrational levels used for the free-bound transi- 
tion, where the resonance position is given in frequency 
detuning from the atomic So — Pi transition. Our setup 
for the PA laser beams could provide enough power for 
detuning smaller than 1 GHz. To obtain enough power 
for detuning larger than 1 GHz, the sideband of the EOM 
was used for the MOT while the carrier is used as the PA 
lasers. We observed the two-color PA signals by measur- 
ing the number of surviving atoms with the absorption 
imaging method with the 1 Sq — X -Pi transition. 



III. EXPERIMENTAL RESULTS 

All the two-color PA spectra observed are shown in 
Fig. [31 where v is the vibrational quantum number count- 
ing down from the dissociation limit and J is the rota- 
tional quantum number. The observed spectral linewidth 
is typically several tens of kHz due to the finite energy 
distribution at the 1 fxK temperature and the power 
broadening. The selection rule for the parity determines 
the observable quantum number J in the ground state 
1 S+ as 0, 2, 4, • • • from the s-wave collision. We observed 
all resonance positions with detuning less than 300 MHz 
from the dissociation limit except for the v = 2 state for 
m Yb; see Fig. 1(b). 

The observed peak positions included the light shift. 



TABLE I: Excited state rovibrational Yb dimer levels used 
for the free-bound PA transition, where v e and J e are the re- 
spective vibrational and rotational quantum numbers; v e is 
numbered from the dissociation limit. The position of the 
excited state level is given as detuning in MHz from the fre- 
quency of the atomic 1 So — 3 Pi transition. For the fermionic 
isotopes the hyperfine level with the largest total atomic an- 
gular momentum for the 3 Pi state defines the dissociation 
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FIG. 4: (color online) Resonance positions /i — fi of the 
v — 1, J = state of 174 Yb2 as a function of 8ls- The 
temperature shift is not compensated. The vertical error bars 
represent the uncertainty estimates for the temperature and 
the center frequency of the resonance. The horizontal error 
bars include the fluctuations of the laser intensities. The data 
are fitted with a linear function represented by a solid line. 
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FIG. 5: (color online) Resonance positions /i — fi of the 
v = 1, J = state of 174 Yb2 as a function of the temperature. 
The light shift is not compensated. The vertical error bars 
represent the uncertainty estimates for the center frequency 
of the resonance. The data are fitted with a linear function 
represented by a solid line. 
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We measured the peak positions with several laser in- 
tensities for Li and L 2 and different detunings A, in or- 
der to compensate the light shift by interpolation or ex- 
trapolation. We found that the bound-bound coupling is 
dominant for the light shift for the Raman spectroscopy. 
Namely, the light shift is expressed as 

^=e(xrhh + i)- (4) 

where j3 is a constant related to the Franck-Condon fac- 
tor of the bound-bound transition, 1\ and I2 are the laser 
intensities of L\ and L 2 , respectively. Typical values of 
the detuning A and laser intensity were about 2 MHz and 
10 — 1000 mW/cm 2 . The contributions of the other ex- 
cited bound states were negligible, since their resonance 
frequencies were typically more than 300 MHz apart, 
which was much larger than A and A + /1 — / 2 . Fig- 
ure[5]shows the light shift 8ls with several intensities for 
L\ and L 2 and detuning A. The error bars include the 
uncertainty due to fluctuations of laser intensity and the 
uncertainty for the temperature of the atom clouds. We 
note that the data are well fitted by a linear function. It 
is also noted that the value of f3 obtained from the fitting 
is consistent with our estimation of the Franck-Condon 
factor of the bound-bound transition, which also ensures 
the validity of our analysis. The resonance position at 
5ls — gives the peak position without the light shift. 
For the Autler-Towncs spectroscopy the light shift was 
basically caused by the laser L\. When the detuning A 
was not zero, however, the laser L2 also contributed to 
the light shift. The sense and magnitude of the detuning 
A could be estimated from the careful investigation of the 
spectral shape (25|. Similarly, the peak positions with 
different laser intensities for L\ and L2 were measured 
and the light shift was compensated by extrapolation of 
these data. 

The observed peak positions also suffered from the 
temperature shift. The temperature shift is assumed to 
be ot^bT 1 , where ar is a constant, fcs is the Boltzmann 
constant and T is the temperature. The factor ax was 
expected to be 3/2 [(|, which was checked experimentally 
for v = 1, J = state of 174 Yb 2 for the temperature range 
from 0.5 /xK to 2 /xK. The resonance positions with dif- 
ferent temperatures are shown in Fig. The data were 
fitted with a linear function represented in a solid line. 
The results are consistent with |fcsT within 9 percent. 
We also compensated the temperature shift for the other 
data using the §/cgT dependence. The final results of the 
experimentally determined binding energies are listed in 
Table HI 

We also measured Eb for the v = 1, J = state of 
174 Yb 2 with the two-color PA using the strongly allowed 
So — Pi transition. The wavelength, lincwidth, and sat- 
uration intensity of the atomic transition were 399 nm, 
29 MHz, and 60 mW/cm 2 , respectively. We used the 
v e \ = 157 level of the excited 4 I]+ molecular state, where 
Vpi is numbered from the dissociation limit Sq + 1 Pi 
fl9l ]. The resonance position of the level was —172.3 



TABLE II: Measured and calculated binding energies Eb for 
homonuclear isotopic pairs, where v and J are the vibrational 
and rotational quantum numbers of the ground state dimer 
level and v is numbered from the dissociation limit. R and AT 
respectively represent the Raman and Autler-Townes spectro- 
scopic method of determining Eb- 



isotope 


V 


J method 
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Difference 










experiment 


theory 


(MHz) 


iYU Yb 


1 





R 


-27.661(23) 


-27.755 


0.094 






2 


R 


-3.651(26) 


-3.683 


0.032 


171 Yb 


1 





AT 


-64.418(40) 


-64.548 


0.130 






2 


AT 


-31.302(50) 


-31.392 


0.090 


172 Yb 


1 





AT 


-123.269(26) 


-123.349 


0.080 






2 


R 


-81.786(19) 


-81.879 


0.093 


173 Yb 


1 





R 


-1.539(74) 


-1.613 


0.074 


174 Yb 


1 





R 


-10.612(38) 


-10.642 


0.030 




1 





AT 


-10.606(17) 


-10.642 


0.036 




2 





R 


-325.607(18) 


-325.607 


0.000 




2 


2 


R 


-268.575(21) 


-268.576 


0.001 


176 yb 


1 





R 


-70.404(11) 


-70.405 


0.001 




1 


2 


R 


-37.142(13) 


-37.118 


-0.024 



GHz from the dissociation limit. Raman spectroscopy 
was employed to measure the peak position. The light 
shift and temperature shift were also compensated. The 
measured Eb of 10.597(59) MHz is in excellent agreement 
with the values in Table [Til which includes the Autler- 
Towncs measurement with a precision of better than 20 
kHz. 



IV. CALCULATION AND DISCUSSION 

The binding energies of the bound states as well 
as scattering lengths of all isotopic combinations are 
determined by the reduced mass and a single Born- 
Oppenheimer potential V(r), as long as small mass- 
dependent adiabatic and non-adiabatic corrections to the 
potential are sufficiently small. The key features of the 
potential that determine the positions of the last few 
bound states are the form of the long range potential 
and a phase associated with the short range potential. 
Consequently, we assume the following simple potential 
form to analyze the data: 

nr) = -^{l-^)-^+B(r)J(J + l), (5) 

where a is a constant, C& is the van der Waals constant 
associated with the dipole-quadrupole interaction, and 
B(r) = h 2 /(2[j,r 2 ) is due to molecular rotation. The first 
term in Eq. ([5]) gives the Lennard- Jones form for the 
potential, for which the short range form can be changed 
by varying a. The C% term is needed to improve the 
quality of the fit to the data. 

By solving the Schrodinger equation numerically for 
the eigenvalues and comparing to the measured binding 
energies, it is possible to determine an optimum set of 



6 



TABLE III: Calculated s-wave scattering lengths in nm for Yb isotopic combinations 
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22.1(7) 
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10.55(11) 
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4.22(10) 


174y b 










5.55(8) 
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potential parameters. The J = 0, v = 2 level of 174 Yb2 
and the J = 0, v = 1 level for 176 Yb2 were fit to de- 
termine Cq and obtain the right number of bound states 
iVi 74 for 174 Yb 2 . Then the J = 2, v = 2 level for 174 Yb 2 
was added to the fit to determine C% and improve the 
determination of C§. The results are Cq = 1931.7 Ehd®, 
C* 8 = 1.93 x 10 5 E h al, a = 9.0109362 a , N 17i = 72, 
where ao « 0.529 nm, Eh f» 4.36 x 10~ 18 J. These param- 
eters then determine without additional adjustment the 
binding energies of the other isotopic combinations shown 
in Table [IT] and Fig [(Jb) . Taking into account a lack of 
precise knowledge about the short range part of the po- 
tential and the magnitude of possible retardation effects 
and neglect of higher order dispersion energy terms, we 
estimate the uncertainties in Cq and Cs to be about 2 per- 
cent and 25 percent respectively, or Cq = 1932(30) Eha^ 
and C 8 = 1.9(5) x 10 5 E h a%. Adding a ~C w /r 10 van der 
Waals term only very slightly improved the quality of 
the fit, and the optimum values of Cq and Cs remain 
within the above stated uncertainties. The Cq value 
is lower than a previous experimental determination of 
2300(250) iS/t^Pil and recent theoretical predictions of 
2291.6 E h a% [26[, 2567.9 E h al H3, and 2062 E h c% [28|. 
The short range form of the potential should be viewed 
as a pscudopotcntial having the right number of bound 
states N but not necessarily giving an accurate shape for 
the potential. The model well depth D e /h = 32469 GHz 
is significantly larger than ab initio values; see Ref. [27j 
and references therein. 

The agreement shown in Table |TT] with a precision less 
than 100 kHz between the experimentally determined Eb 
and most of the calculated Eb values is quite impressive 
for such a simple model potential. The use of a single 
mass-independent potential with the appropriate reduced 
mass is thus seen as an excellent approximation for cal- 
culating the isotopic variation in binding energies. The 
failure to obtain a fit to the data within experimental er- 
ror in all cases could be indicative of the failure of mass 
scaling, although it may only be due to the limitations 
of the simple form we assumed to represent the potential 
over its whole range. Adding a small mass-dependent 
correction on the order of 1 GHz to the well depth of the 
potential allows us to fit the binding energies for each iso- 
tope almost within experimental uncertainties. We sec 
indications that the potential is deeper for heavier iso- 



topes. Additional work is needed to see if the Yb system 
can be used to make quantitative tests of the accuracy 
and limitations of mass scaling. 

When we added terms to account for relativistic re- 
tardation effects [lU, the \ 2 of the overall fit did not 
improve, although the energies of the least bound levels 
improved slightly. These relativistic terms take on the 
form a 2 W4/r i in the physically interesting region r -C A, 
where a is the fine structure constant and X ~ 1100 ao is 
a characteristic length associated with a mean electronic 
excitation energy. The lead term in the retarded van der 
Waals interactions switches to 1/r 7 behavior in the very 
long range region r»A. Since the outer turning point 
of the least bound level is still only around 150 ao, the 
o^H^/r 4 form is the appropriate form to use in looking 
for the magnitude of the effect. When we added a C^/r 4 
term to the potential, we find that we get a slightly bet- 
ter fit for the most weakly bound levels, but not a better 
overall fit, if we take C4 = 1.5 x 10~ 3 Ehd^, similar in 
magnitude to the value 0.6 x 10~ 3 E^a^ estimated in the 
way proposed in Ref. p9l | using the dipole polarizability 
and Cq coefficient for Yb reported in Ref. [26|. Our re- 
ported model parameters and scattering lengths include 
the uncertainties associated with the lack of knowledge 
of the retardation corrections for this system. Conse- 
quently, there is a need for a theoretical evaluation of 
these corrections. 

Given the ground state potential energy curve with the 
form of Eq. ([5]) and the parameters given above, we can 
calculate the s-wave scattering lengths for all isotopic 
combinations by numerically solving the Schrodingcr 
equation with the appropriate reduced mass. These are 
shown in Table Hill and Fig. [6]^a). The uncertainties re- 
flect the uncertainties in the model parameters and the 
need to include retardation corrections to the potential. 
The value for 174 Yb is in excellent agreement with the 
value 5.53(11) reported in Ref. [lj]. The value for 173 Yb 
is slightly larger than the value estimated from the ther- 
malization experiments in Ref. [bsj . However, the s-wave 
elastic cross section for 173 Yb decreases by about 10 per- 
cent between E = and the 6 /xK temperature of the 
experiment. The experimental value could also be larger 
if the unknown distribution of spin populations differed 
from the assumption of uniformity. 

Since the three-parameter model potential Eq. (JSJ) is 
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FIG. 6: (color online) Calculated scattering lengths (a) and 
binding energies Eh (b) versus twice the reduced mass, 2/i in 
the atomic mass unit, using the 3-parameter potential energy 
model in Eq. ([5} with the parameters given in the text. The 
vertical lines show the masses for the 7 like-atom pairs. The 
solid and dashed lines in (b) show the J — and 2 eigenvalues 
respectively. The measured values with vertical error bars are 
shown for the levels for which they have been measured. The 
measured 174 Yb scattering length of Ref. is also shown 
in (a). The horizontal dashed line in (a) shows the van der 
Waals length a. 



assumed to be common for all possible isotopic combina- 
tions, the simple analytical formula Eq. (fTJ) also holds. 
The predictions of this formula are completely indistin- 
guishable from those of the numerical calculation on the 
scale of Fig. [S] (a). The actual difference between scatter- 
ing lengths calculated exactly from the model potential 
and those calculated from the analytical formula in Eq. 
([1]) are below 0.03 nm for all like isotope cases except the 
one with the largest scattering length magnitude, 172 Yb, 
for which the difference is 0.35 nm. A similar statement 
applies to the mixed-isotope cases, for which most differ- 
ences are below 0.03 nm. 

As the reduced mass varies from 168/2 to 176/2, the 
scattering length varies through a complete cycle from 
—oo to +oo. In fact, we have a rich variety of the scatter- 
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FIG. 7: (color online) Calculated cross section a(E) versus 
collision energy Ejks for two atoms of the same isotope. The 
label for each curve shows the isotopic mass number. Even 
mass number corresponds to identical bosons and odd mass 
number corresponds to two fermions with different spin pro- 
jections. The dashed lines show the s-wave contribution to 
the total cross section. The solid lines show the contribution 
from s- and d- waves for the like boson cases, and the con- 
tribution from s-, p-, and d-waves for the fermion case with 
different spin components. The E — > cross sections are 87ra 2 
and 47ra 2 for the respective boson and fermion cases. 



ing lengths from large negative values for 172 Yb - 172 Yb, 

171 Yb _173 Ybj 170 Yb _174 Yb; and 168 Yb _176 Yb; &nd &1 _ 

most zero for 171 Yb - 171 Yb, 170 Yb - 172 Yb, and 168 Yb 
- 174 Yb, and to large positive for 172 Yb - 173 Yb and 171 Yb 
_i74 Yb Thus, t be scattering len gth can be widely tuned 
by varying isotopic composition [30(. 

It should be noted that the observed behavior of evap- 
orative cooling and sympathetic cooling arc consistent 
with these scattering lengths. The efficient evaporative 
cooling for 170 Yb, 173 Yb, and 174 Yb is consistent with 
the large scattering lengths: 3.4 nm for 170 Yb, 11 nm 
for 173 Yb, and 5.6 nm for 174 Yb. Inefficient evaporative 
cooling for 171 Yb and 176 Yb is also consistent with the 
small scattering lengths —0.1 nm for 171 Yb and —1.3 nm 
for 176 Yb. On the other hand, we have successfully used 
sympathetic cooling to cool 171 Yb and 176 Yb with 174 Yb, 
for which the respective mixed-isotope scattering lengths 
are 22 nm and 2.9 nm. The extremely large value of —32 
nm for 172 Yb would explain very rapid atom decay ob- 
served for this system, which could be due to three-body 
recombination. 

Using the three-parameter model potential, we can also 
calculate the collisional properties of other partial waves 
at non-zero collision energies. Figure [7] shows the energy- 
dependent cross section for the collision of like isotopic 
species. The cross section can be resonantly enhanced by 
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a shape resonance, caused by the existence of quasibound 
rovibrational levels supported by the centrifugal barrier. 
A d-wave shape resonance exists for 174 Yb as already 
pointed out in Refs. 0, [H|. Our model predicts a broad 
peak in the collision cross section near E/h = 4.5 MHz, 
or E/ks = 220 (jK, which is off-scale in Fig. [7] We found 
that a low energy p-wave shape resonance also exists for 
173 Yb, giving rise to a peak in the cross section near 
E/k B = 48 fiK. 

The s-wave contribution to the cross section becomes 
zero when the s-wave collisional phase shift is zero. Such 
a zero as a function of energy results in a Ramsauer- 
Townsend minimum in the cross section. This minimum 
is especially evident when the collision energy is so small 
that other partial waves make negligible contributions to 
the cross section, This effect occurs at very low collision 
energy when the scattering length has a small negative 
value, as for 171 Yb and for 176 Yb. Figure shows the 
calculated cross section minima near 2 and 25 /j,K for 
these respective species. This effect explains why evapo- 
rative cooling is found to be inefficient for these isotopes. 

V. CONCLUDING REMARKS 

In conclusion, we report the accurate determination 
of the binding energies of the twelve least bound states 
in the ground molecular potentials for six Yb isotopes 
and the accurate determination of the s-wave scattering 
lengths for all possible combination of the isotopes based 
on a simple three-parameter model potential. The model 



parameters are based on fitting a very limited set of the 
experimental data, making use of only three binding en- 
ergies: two bound-states of 174 Yb2 and one bound-state 
of 176 Yb2- Fitting this limited set of data allows us to 
reproduce the other nine binding energies with an accu- 
racy of about 100 kHz and to determine the 28 scatter- 
ing lengths for the different isotope combinations with 
an accuracy of about a few percent for most cases. In 
addition, we can calculate the energy-dependent elastic 
scattering due to other partial waves such as p- and d- 
waves. These results provide an important foundation 
for future research with Yb atoms on such topics as the 
efficiency of evaporative cooling, the stability of quantum 
gases and their mixtures, and clock shifts. 
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